Return to MUSA 801
Projects Page
This project was completed for the MUSA/Smart Cities Practicum course
(MUSA 801) instructed by Michael Fichman and Matthew Harris. We are
grateful to our instructors for their continued support and feedback. We
would like to give special thanks to KC Filippino and Ben McFarlane from
Hampton Roads Planning District Commission, and Dexter Locke from the
United States Forest Service for providing data, insight, and support
throughout the semester. This project would not have been possible
without them.
1.Introduction
1.1 Abstract
This project aims to develop a precision forecast model for land
cover change at the Chesapeake Watershed, focusing on three
representative counties: Isle of Wight, James City, and Portsmouth. By
leveraging high-resolution longitudinal land use and land cover data
provided by the Chesapeake Conservancy, the model will predict land
cover conversions from pervious to impervious surfaces. This forecast
will enable land use and environmental planners to visualize and assess
potential impacts on flood risk, heat stress, and heat island risks in
the region. The model will be generalizable to the county level,
incorporating only widely available inputs, thus allowing any
municipality within the Chesapeake basin to replicate the analysis. This
proof-of-concept project will demonstrate the utility of precision
conservation in climate adaptation and mitigation planning and provide a
valuable tool for planners and policymakers across the region.

1.2 Motivation & Use Case
To build resilient communities, the HRPDC has set up a green
infrastructure plan. This plan aims to to identify and prioritize a
network of valuable conservation lands in order to achieve multiple
benefits, such as habitat protection, drinking water supply protection,
stormwater management and recreational opportunities. A new component of
this plan is to build a model for the potential future growth and
identify which areas of the green infrastructure network are most at
risk for development.
2. Exploratory Analysis
2.1 Understanding landcover data
The Chesapeake Conservancy supplies us with high-resolution landcover
data, which is essential for our precision conservation efforts. This
data is a vast raster dataset with an impressive 1-meter accuracy,
offering 900 times more detail than the commonly used 30-meter
resolution National Land Cover Dataset. This level of detail is critical
in capturing subtle changes in land use and land cover. Within the
landcover classification, pervious surfaces include categories such as
tree canopy, shrub, and wetlands, which allow water to infiltrate the
ground, while impervious surfaces encompass categories like roads and
structures that prevent water infiltration, leading to increased runoff
and potential flooding issues. This detailed classification enables us
to better understand and predict land cover changes, particularly the
conversion from pervious to impervious surfaces.

Data
Source
2.2 Data Cleaning and Wrangling
2.3
dependent: pervious to impervious surface during the period
2014-201
independent:
lccover change 能change
1. `originallc` originallandcover
2. `lcp` permeable/impermeable:
< 6 - 0
>= 6 - 1
3. `lcchange` whether landcover has changed: 1
4. `lc` 3*3的permeable rate:
0-1(0.1,0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8)
5. `popchange` per m^2
6. `pctwhitechange`
7. `unitchange` per m^2
8. `medhhincchange` per m^2
9. `road` road 1/0
10. `water` water
11. `canopy` tree canopy:
tree canopy - 1
tree canopy over structure -0.5
tree canopy over other -0.5
tree canopy over road - 0.5
other - 0
12. `perm` permeable:
2 -shrub - 0.3
4 - herb - 0.3
5 - wetand - 0.4
other -0
13. `barren` barren:
6 - barren - 1
not barren - 0
14. `struct` structure:
7 - structure - 1
not structure - 0
15. `other` other:
8 - other - 1
not other - 0
16. `terrain`dem
17. `slope`slope
18. `area` m^2 per block group
19. `geoid`-- blockgroup id
20. `soil type`-
模型encode:categorial:尤其是
possible features per 10*10 :
Road density distance to rivers
impervious / pervious landcover type population density population
change per 4 year median hh income median hh income change per 4 year
percentage of white population percentage of white population change per
4 year percentage of changing from impervious to pervious/from pervious
to impervious in 4 years predicted landcover type:
possible features per block group :
percentage of impervious / pervious change: population population
change per 4 year median hh income median hh income change per 4 year
percentage of white population percentage of white population change per
4 year percentage of changing from impervious to pervious/from pervious
to impervious in 4 years predicted landcover type composition:
possible feautres per county:
percentage of impervious / pervious change: population population
change per 4 year median hh income median hh income change per 4 year
percentage of white population percentage of white population change per
4 year percentage of changing from impervious to pervious/from pervious
to impervious in 4 years predicted landcover type composition:
development speed defined by: speed of population change and landcover
change speed Agricultural land use, vegetation cover and water loss:
LS0tDQp0aXRsZTogIlByZWNpc2lvbiBGb3JlY2FzdHMgb2YgTGFuZCBDb3ZlciBDaGFuZ2UiDQpzdWJ0aXRsZTogIkNoZXNhcGVha2UgV2F0ZXJzaGVkIg0KYXV0aG9yOiAiWXVld2VuIERhaSwgU2h1amluZyBZaSwgWGluZ2UgWmhhbmciDQpkYXRlOiAiMjAyMy0wNC0yNSINCm91dHB1dDogDQogIGh0bWxfZG9jdW1lbnQ6DQogICAgdG9jOiB0cnVlDQogICAgdG9jX2Zsb2F0OiB0cnVlDQogICAgY29kZV9mb2xkaW5nOiAiaGlkZSINCiAgICBjb2RlX2Rvd25sb2FkOiB0cnVlDQotLS0NCg0KYGBge3Igc2V0dXAsIGluY2x1ZGU9RkFMU0V9DQprbml0cjo6b3B0c19jaHVuayRzZXQoZWNobyA9IFRSVUUpDQpgYGANCg0KW1JldHVybiB0byBNVVNBIDgwMSBQcm9qZWN0cyBQYWdlXShodHRwczovL3Blbm5tdXNhLmdpdGh1Yi5pby9NVVNBXzgwMS5pby8pDQoNClRoaXMgcHJvamVjdCB3YXMgY29tcGxldGVkIGZvciB0aGUgTVVTQS9TbWFydCBDaXRpZXMgUHJhY3RpY3VtIGNvdXJzZSAoTVVTQSA4MDEpIGluc3RydWN0ZWQgYnkgTWljaGFlbCBGaWNobWFuIGFuZCBNYXR0aGV3IEhhcnJpcy4gV2UgYXJlIGdyYXRlZnVsIHRvIG91ciBpbnN0cnVjdG9ycyBmb3IgdGhlaXIgY29udGludWVkIHN1cHBvcnQgYW5kIGZlZWRiYWNrLiBXZSB3b3VsZCBsaWtlIHRvIGdpdmUgc3BlY2lhbCB0aGFua3MgdG8gS0MgRmlsaXBwaW5vIGFuZCBCZW4gTWNGYXJsYW5lIGZyb20gSGFtcHRvbiBSb2FkcyBQbGFubmluZyBEaXN0cmljdCBDb21taXNzaW9uLCBhbmQgRGV4dGVyIExvY2tlIGZyb20gdGhlIFVuaXRlZCBTdGF0ZXMgRm9yZXN0IFNlcnZpY2UgZm9yIHByb3ZpZGluZyBkYXRhLCBpbnNpZ2h0LCBhbmQgc3VwcG9ydCB0aHJvdWdob3V0IHRoZSBzZW1lc3Rlci4gVGhpcyBwcm9qZWN0IHdvdWxkIG5vdCBoYXZlIGJlZW4gcG9zc2libGUgd2l0aG91dCB0aGVtLg0KDQojIyAxLkludHJvZHVjdGlvbg0KIyMjIDEuMSBBYnN0cmFjdA0KVGhpcyBwcm9qZWN0IGFpbXMgdG8gZGV2ZWxvcCBhIHByZWNpc2lvbiBmb3JlY2FzdCBtb2RlbCBmb3IgbGFuZCBjb3ZlciBjaGFuZ2UgYXQgdGhlIENoZXNhcGVha2UgV2F0ZXJzaGVkLCBmb2N1c2luZyBvbiB0aHJlZSByZXByZXNlbnRhdGl2ZSBjb3VudGllczogSXNsZSBvZiBXaWdodCwgSmFtZXMgQ2l0eSwgYW5kIFBvcnRzbW91dGguIEJ5IGxldmVyYWdpbmcgaGlnaC1yZXNvbHV0aW9uIGxvbmdpdHVkaW5hbCBsYW5kIHVzZSBhbmQgbGFuZCBjb3ZlciBkYXRhIHByb3ZpZGVkIGJ5IHRoZSBDaGVzYXBlYWtlIENvbnNlcnZhbmN5LCB0aGUgbW9kZWwgd2lsbCBwcmVkaWN0IGxhbmQgY292ZXIgY29udmVyc2lvbnMgZnJvbSBwZXJ2aW91cyB0byBpbXBlcnZpb3VzIHN1cmZhY2VzLiBUaGlzIGZvcmVjYXN0IHdpbGwgZW5hYmxlIGxhbmQgdXNlIGFuZCBlbnZpcm9ubWVudGFsIHBsYW5uZXJzIHRvIHZpc3VhbGl6ZSBhbmQgYXNzZXNzIHBvdGVudGlhbCBpbXBhY3RzIG9uIGZsb29kIHJpc2ssIGhlYXQgc3RyZXNzLCBhbmQgaGVhdCBpc2xhbmQgcmlza3MgaW4gdGhlIHJlZ2lvbi4gVGhlIG1vZGVsIHdpbGwgYmUgZ2VuZXJhbGl6YWJsZSB0byB0aGUgY291bnR5IGxldmVsLCBpbmNvcnBvcmF0aW5nIG9ubHkgd2lkZWx5IGF2YWlsYWJsZSBpbnB1dHMsIHRodXMgYWxsb3dpbmcgYW55IG11bmljaXBhbGl0eSB3aXRoaW4gdGhlIENoZXNhcGVha2UgYmFzaW4gdG8gcmVwbGljYXRlIHRoZSBhbmFseXNpcy4gVGhpcyBwcm9vZi1vZi1jb25jZXB0IHByb2plY3Qgd2lsbCBkZW1vbnN0cmF0ZSB0aGUgdXRpbGl0eSBvZiBwcmVjaXNpb24gY29uc2VydmF0aW9uIGluIGNsaW1hdGUgYWRhcHRhdGlvbiBhbmQgbWl0aWdhdGlvbiBwbGFubmluZyBhbmQgcHJvdmlkZSBhIHZhbHVhYmxlIHRvb2wgZm9yIHBsYW5uZXJzIGFuZCBwb2xpY3ltYWtlcnMgYWNyb3NzIHRoZSByZWdpb24uDQoNCg0KIVtdKGltYWdlcy9IUlBEQy5qcGcpDQoNCg0KIyMjIDEuMiBNb3RpdmF0aW9uICYgVXNlIENhc2UNCg0KVG8gYnVpbGQgcmVzaWxpZW50IGNvbW11bml0aWVzLCB0aGUgSFJQREMgaGFzIHNldCB1cCBhIGdyZWVuIGluZnJhc3RydWN0dXJlIHBsYW4uIFRoaXMgcGxhbiBhaW1zIHRvIHRvIGlkZW50aWZ5IGFuZCBwcmlvcml0aXplIGEgbmV0d29yayBvZiB2YWx1YWJsZSBjb25zZXJ2YXRpb24gbGFuZHMgaW4gb3JkZXIgdG8gYWNoaWV2ZSBtdWx0aXBsZSBiZW5lZml0cywgc3VjaCBhcyBoYWJpdGF0IHByb3RlY3Rpb24sIGRyaW5raW5nIHdhdGVyIHN1cHBseSBwcm90ZWN0aW9uLCBzdG9ybXdhdGVyIG1hbmFnZW1lbnQgYW5kIHJlY3JlYXRpb25hbCBvcHBvcnR1bml0aWVzLiBBIG5ldyBjb21wb25lbnQgb2YgdGhpcyBwbGFuIGlzIHRvIGJ1aWxkIGEgbW9kZWwgZm9yIHRoZSBwb3RlbnRpYWwgZnV0dXJlIGdyb3d0aCBhbmQgaWRlbnRpZnkgd2hpY2ggYXJlYXMgb2YgdGhlIGdyZWVuIGluZnJhc3RydWN0dXJlIG5ldHdvcmsgYXJlIG1vc3QgYXQgcmlzayBmb3IgZGV2ZWxvcG1lbnQuDQoNCiMjIDIuIEV4cGxvcmF0b3J5IEFuYWx5c2lzDQojIyMgMi4xIFVuZGVyc3RhbmRpbmcgbGFuZGNvdmVyIGRhdGENClRoZSBDaGVzYXBlYWtlIENvbnNlcnZhbmN5IHN1cHBsaWVzIHVzIHdpdGggaGlnaC1yZXNvbHV0aW9uIGxhbmRjb3ZlciBkYXRhLCB3aGljaCBpcyBlc3NlbnRpYWwgZm9yIG91ciBwcmVjaXNpb24gY29uc2VydmF0aW9uIGVmZm9ydHMuIFRoaXMgZGF0YSBpcyBhIHZhc3QgcmFzdGVyIGRhdGFzZXQgd2l0aCBhbiBpbXByZXNzaXZlIDEtbWV0ZXIgYWNjdXJhY3ksIG9mZmVyaW5nIDkwMCB0aW1lcyBtb3JlIGRldGFpbCB0aGFuIHRoZSBjb21tb25seSB1c2VkIDMwLW1ldGVyIHJlc29sdXRpb24gTmF0aW9uYWwgTGFuZCBDb3ZlciBEYXRhc2V0LiBUaGlzIGxldmVsIG9mIGRldGFpbCBpcyBjcml0aWNhbCBpbiBjYXB0dXJpbmcgc3VidGxlIGNoYW5nZXMgaW4gbGFuZCB1c2UgYW5kIGxhbmQgY292ZXIuIFdpdGhpbiB0aGUgbGFuZGNvdmVyIGNsYXNzaWZpY2F0aW9uLCBwZXJ2aW91cyBzdXJmYWNlcyBpbmNsdWRlIGNhdGVnb3JpZXMgc3VjaCBhcyB0cmVlIGNhbm9weSwgc2hydWIsIGFuZCB3ZXRsYW5kcywgd2hpY2ggYWxsb3cgd2F0ZXIgdG8gaW5maWx0cmF0ZSB0aGUgZ3JvdW5kLCB3aGlsZSBpbXBlcnZpb3VzIHN1cmZhY2VzIGVuY29tcGFzcyBjYXRlZ29yaWVzIGxpa2Ugcm9hZHMgYW5kIHN0cnVjdHVyZXMgdGhhdCBwcmV2ZW50IHdhdGVyIGluZmlsdHJhdGlvbiwgbGVhZGluZyB0byBpbmNyZWFzZWQgcnVub2ZmIGFuZCBwb3RlbnRpYWwgZmxvb2RpbmcgaXNzdWVzLiBUaGlzIGRldGFpbGVkIGNsYXNzaWZpY2F0aW9uIGVuYWJsZXMgdXMgdG8gYmV0dGVyIHVuZGVyc3RhbmQgYW5kIHByZWRpY3QgbGFuZCBjb3ZlciBjaGFuZ2VzLCBwYXJ0aWN1bGFybHkgdGhlIGNvbnZlcnNpb24gZnJvbSBwZXJ2aW91cyB0byBpbXBlcnZpb3VzIHN1cmZhY2VzLg0KDQohW10oaW1hZ2VzL0RhdGEucG5nKQ0KDQpbRGF0YSBTb3VyY2VdKGh0dHBzOi8vd3d3LmNoZXNhcGVha2Vjb25zZXJ2YW5jeS5vcmcvY29uc2VydmF0aW9uLWlubm92YXRpb24tY2VudGVyL2hpZ2gtcmVzb2x1dGlvbi1kYXRhL2x1bGMtZGF0YS1wcm9qZWN0LTIwMjIvKQ0KDQoNCiMjIyAyLjIgRGF0YSBDbGVhbmluZyBhbmQgV3JhbmdsaW5nDQoNCiMjIyAyLjMgDQoNCiMjIyMgZGVwZW5kZW50OiBwZXJ2aW91cyB0byBpbXBlcnZpb3VzIHN1cmZhY2UgZHVyaW5nIHRoZSBwZXJpb2QgMjAxNC0yMDENCiMjIyMgaW5kZXBlbmRlbnQ6DQogICAgbGNjb3ZlciBjaGFuZ2Ug6IO9Y2hhbmdlDQogICAgMS4gYG9yaWdpbmFsbGNgIG9yaWdpbmFsbGFuZGNvdmVyDQogICAgDQogICAgMi4gYGxjcGAgcGVybWVhYmxlL2ltcGVybWVhYmxlOg0KICAgIDwgNiAtIDANCiAgICA+PSA2IC0gMQ0KICAgIA0KICAgIDMuIGBsY2NoYW5nZWAgd2hldGhlciBsYW5kY292ZXIgaGFzIGNoYW5nZWQ6IDENCiAgICANCiAgICA0LiBgbGNgIDMqM+eahHBlcm1lYWJsZSByYXRlOg0KICAgIDAtMe+8iDAuMe+8jDAuMu+8jCAwLjPvvIwgMC4077yMIDAuNe+8jCAwLjbvvIwgMC4377yMIDAuOO+8iQ0KICAgIA0KICAgIDUuIGBwb3BjaGFuZ2VgIHBlciBtXjINCiAgICANCiAgICA2LiBgcGN0d2hpdGVjaGFuZ2VgDQogICAgDQogICAgNy4gYHVuaXRjaGFuZ2VgIHBlciBtXjINCiAgICANCiAgICA4LiBgbWVkaGhpbmNjaGFuZ2VgIHBlciBtXjINCiAgICANCiAgICA5LiBgcm9hZGAgcm9hZCAxLzANCiAgICANCiAgICAxMC4gYHdhdGVyYCB3YXRlcg0KICAgIA0KICAgIDExLiBgY2Fub3B5YCB0cmVlIGNhbm9wee+8mg0KICAgIHRyZWUgY2Fub3B5IC0gMQ0KICAgIHRyZWUgY2Fub3B5IG92ZXIgc3RydWN0dXJlIC0wLjUNCiAgICB0cmVlIGNhbm9weSBvdmVyIG90aGVyIC0wLjUNCiAgICB0cmVlIGNhbm9weSBvdmVyIHJvYWQgLSAwLjUNCiAgICBvdGhlciAtIDANCiAgICANCiAgICAxMi4gYHBlcm1gIHBlcm1lYWJsZToNCiAgICAyIC1zaHJ1YiAtIDAuMw0KICAgIDQgLSBoZXJiIC0gMC4zDQogICAgNSAtIHdldGFuZCAtIDAuNA0KICAgIG90aGVyIC0wDQogICAgDQogICAgMTMuIGBiYXJyZW5gIGJhcnJlbjoNCiAgICA2IC0gYmFycmVuIC0gMQ0KICAgIG5vdCBiYXJyZW4gLSAwDQogICAgDQogICAgMTQuIGBzdHJ1Y3RgIHN0cnVjdHVyZToNCiAgICA3IC0gc3RydWN0dXJlIC0gMQ0KICAgIG5vdCBzdHJ1Y3R1cmUgLSAwDQogICAgDQogICAgMTUuIGBvdGhlcmAgb3RoZXI6DQogICAgOCAtIG90aGVyIC0gMQ0KICAgIG5vdCBvdGhlciAtIDANCiAgICANCiAgICAxNi4gYHRlcnJhaW5gZGVtDQogICAgDQogICAgMTcuIGBzbG9wZWBzbG9wZQ0KICAgIA0KICAgIDE4LiBgYXJlYWAgbV4yIHBlciBibG9jayBncm91cA0KICAgIA0KICAgIDE5LiBgZ2VvaWRgLS0gYmxvY2tncm91cCBpZA0KICAgIA0KICAgIDIwLiBgc29pbCB0eXBlYC0NCg0KDQrmqKHlnotlbmNvZGXvvJpjYXRlZ29yaWFs77ya5bCk5YW25pivDQoNCnBvc3NpYmxlIGZlYXR1cmVzIHBlciAxMCoxMCA6DQoNClJvYWQgZGVuc2l0eQ0KKipkaXN0YW5jZSB0byByaXZlcnMqKg0KDQppbXBlcnZpb3VzIC8gcGVydmlvdXMNCmxhbmRjb3ZlciB0eXBlDQpwb3B1bGF0aW9uIGRlbnNpdHkNCnBvcHVsYXRpb24gY2hhbmdlIHBlciA0IHllYXINCm1lZGlhbiBoaCBpbmNvbWUNCm1lZGlhbiBoaCBpbmNvbWUgY2hhbmdlIHBlciA0IHllYXINCnBlcmNlbnRhZ2Ugb2Ygd2hpdGUgcG9wdWxhdGlvbg0KcGVyY2VudGFnZSBvZiB3aGl0ZSBwb3B1bGF0aW9uIGNoYW5nZSBwZXIgNCB5ZWFyDQpwZXJjZW50YWdlIG9mIGNoYW5naW5nIGZyb20gaW1wZXJ2aW91cyB0byBwZXJ2aW91cy9mcm9tIHBlcnZpb3VzIHRvIGltcGVydmlvdXMgaW4gNCB5ZWFycw0KcHJlZGljdGVkIGxhbmRjb3ZlciB0eXBlOg0KDQpwb3NzaWJsZSBmZWF0dXJlcyBwZXIgYmxvY2sgZ3JvdXAgOg0KDQpwZXJjZW50YWdlIG9mIGltcGVydmlvdXMgLyBwZXJ2aW91cyBjaGFuZ2U6DQpwb3B1bGF0aW9uDQpwb3B1bGF0aW9uIGNoYW5nZSBwZXIgNCB5ZWFyDQptZWRpYW4gaGggaW5jb21lDQptZWRpYW4gaGggaW5jb21lIGNoYW5nZSBwZXIgNCB5ZWFyDQpwZXJjZW50YWdlIG9mIHdoaXRlIHBvcHVsYXRpb24NCnBlcmNlbnRhZ2Ugb2Ygd2hpdGUgcG9wdWxhdGlvbiBjaGFuZ2UgcGVyIDQgeWVhcg0KcGVyY2VudGFnZSBvZiBjaGFuZ2luZyBmcm9tIGltcGVydmlvdXMgdG8gcGVydmlvdXMvZnJvbSBwZXJ2aW91cyB0byBpbXBlcnZpb3VzIGluIDQgeWVhcnMNCnByZWRpY3RlZCBsYW5kY292ZXIgdHlwZSBjb21wb3NpdGlvbjoNCg0KcG9zc2libGUgZmVhdXRyZXMgcGVyIGNvdW50eToNCg0KcGVyY2VudGFnZSBvZiBpbXBlcnZpb3VzIC8gcGVydmlvdXMgY2hhbmdlOg0KcG9wdWxhdGlvbg0KcG9wdWxhdGlvbiBjaGFuZ2UgcGVyIDQgeWVhcg0KbWVkaWFuIGhoIGluY29tZQ0KbWVkaWFuIGhoIGluY29tZSBjaGFuZ2UgcGVyIDQgeWVhcg0KcGVyY2VudGFnZSBvZiB3aGl0ZSBwb3B1bGF0aW9uDQpwZXJjZW50YWdlIG9mIHdoaXRlIHBvcHVsYXRpb24gY2hhbmdlIHBlciA0IHllYXINCnBlcmNlbnRhZ2Ugb2YgY2hhbmdpbmcgZnJvbSBpbXBlcnZpb3VzIHRvIHBlcnZpb3VzL2Zyb20gcGVydmlvdXMgdG8gaW1wZXJ2aW91cyBpbiA0IHllYXJzDQpwcmVkaWN0ZWQgbGFuZGNvdmVyIHR5cGUgY29tcG9zaXRpb246DQpkZXZlbG9wbWVudCBzcGVlZCBkZWZpbmVkIGJ5OiBzcGVlZCBvZiBwb3B1bGF0aW9uIGNoYW5nZSBhbmQgbGFuZGNvdmVyIGNoYW5nZSBzcGVlZA0KQWdyaWN1bHR1cmFsIGxhbmQgdXNlLCB2ZWdldGF0aW9uIGNvdmVyIGFuZCB3YXRlciBsb3NzOg0KDQoNCg0K